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Abstract 

A toy top is defined as a rotationally symmetric body moving in a constant gravita- 
tional field while one point on the symmetry axis is constrained to stay in a horizontal 
plane. It is an integrable system similar to the Lagrange top. Euler-Poisson equa- 
tions are derived. Following Felix Klein, the special unitary group SU(2) is used as 
configuration space and the solution is given in terms of hyperelliptic integrals. The 
curve traced by the point moving in the horizontal plane is analyzed, and a qualitative 
classification is achieved. The cases in which the hyperelliptic integrals degenerate to 
elliptic ones are found and the corresponding solutions are given in terms of Weier- 
strass elliptic functions. 

1 Introduction 

The three famous integrable cases of rigid body motion, the tops of Euler, Lagrange and 
Kowalewsky, have been paramount examples in the theory of integrable systems. The 
modern algebra-geometric approach, using Lax pairs with a spectral parameter 
has been applied to all three. It is surprising that the following system appears only 
sporadically in the literature: A rotationally symmetric rigid body moving in a homoge- 
neous gravitational field with one point on its axis not fixed, but constrained to move in 
a horizontal plane. Following F. Klein p. 58], we call such a system a toy top. 

In many ways, the toy top is similar to Lagrange's top. It is completely integrable 
due to the same kind of rotational symmetry. The solution leads to hyperelliptic integrals 
instead of elliptic ones, but their analytic properties are similar. It seems that Poisson 
is the first to solve the system Q, using Euler angles. It is treated similarly by E. T. 
Whittaker and F. Klein Q. Later, Klein discovered that, as in the case of Lagrange's 
top, simpler solutions are obtained when SU(2) is used as configuration space 0) Hi- 

In this paper we close some gaps in the classical treatment. Following A. I. Bobenko's 
and Yu. B. Suris' treatment of the Lagrange top [p, Q, new equations of motion are 
derived within the framework of Lagrangian mechanics on Lie groups. 
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As the toy top moves, its tip traces a curve on the supporting plane. Formulas for this 
curve are derived and qualitatively different cases are classified. 

For certain values of the first integrals, the hyperelliptic integrals appearing in the 
solution degenerate to elliptic integrals. These cases are classified and the corresponding 
solutions are given in terms of Weierstrass elliptic functions. 



2 Preliminaries 

After it is shown how the group of rotations can be considered the configuration space 
of the toy top, the alternative use of SU(2) is discussed. Finally, Mobius transformations 
that rotate the Riemann sphere of numbers are discussed for use in chapter ^ 



2.1 The group SU(2) as configuration space 

The configuration space of the toy top is the direct product of the rotation group with the 
group of translations in the plane. But since both the gravitational and the resistive force 
of the plane are vertical, the horizontal component of the velocity of the top's center of 
mass is constant. After changing to a suitable moving coordinate system if necessary, the 
center of mass moves only vertically. This reduces the configuration space to the rotation 
group: The top can be brought into any admissible position by a rotation around its center 
of mass, followed by a vertical translation. But the latter is determined by the former due 
to the constraint. 

Instead of the special orthogonal group S0(3) or Euler angles, the special unitary group 
SU(2) will be used to describe configurations of the toy top. As F. Klein discovered, this 
leads to simpler formulas. The special unitary group of two dimensions is the matrix group 



SU(2) 



7 6 



a6 — (3j = 1, S = a, 7 = —(3 



Its Lie algebra consists of the skew hermitian matrices with trace zero: 



su(2) 



1x5 



^X2 + IXi 

Identify M"^ with su(2) via 



-X2 + ixi 
-1x3 



G C 



2x2 



Xi,X2,X3 G 



{xi,X2,X3) 



1 / iX3 -X2 + ixi 

2 1x2 + ixi —ix3 



This corresponds to the choice of basis 



1 /O i 



'' = 2U 



1 /O 



2 Vl 



1 fi 



"=^=2 VO 



The factor 1/2 is introduced so that the cross product in corresponds to the Lie bracket 
in su(2): 



X xy = [x,y]. 
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The induced scalar product on su(2) is 
{x,y) = -2Tr(xy). 

The adjoint action of SU(2) on its Lie algebra is orthogonal and orientation preserving. 
The homomorphism 

Ad : SU(2) ^ S0(3) 

$ ^ (Ad# : X ^ (^x^~^) 

is 2 to 1 with Kernel {±1}- 

Consider two orthonormal coordinate systems whose origin is the center of mass of the 
toy top. The first, called the fixed frame, has axes whose directions are fixed. Assume that 
the third axis points upward. The second coordinate system, called the moving frame, 
moves with the toy top. Assume its third axis points along the top's symmetry axis away 
from the top's tip. The momentary position of the top is given by a matrix $ G SU(2) 
such that \i X = {Xi,X2,X^) is the coordinate vector of some point in the body frame 
then its coordinates x = {xi,X2,x^) with respect to the fixed frame are given by 

1 ( ix3 -X2 + ixi\ 1^/^ iX^ -X2 + iXi\ I 



, ^ ^ ' ' ' $-1. (2.1) 

2 yx2 + ixi -ix3 J 2 V^2 + iXi -iX^ J ^ ' 

In the context of rigid body mechanics, the entries of an SU(2) matrix are called Cayley- 
Klein parameters. 

Suppose that a curve in SU(2) describing the motion of the top is given by 
/a(t) Pit) 

A moving point whose coordinate vector is X(t) in the moving frame has coordinates x{t) = 
^{t)X{t)^^^{t) in the fixed frame. Taking the derivative, one obtains x' = ^X'^~^ + [uj, x] 
and X' = ^^^x'^ - [0.,X], where ui = $'$-1 and Q = <I>~^<I>'. Since ui = they 
represent the same vector, the angular velocity vector, in the fixed and moving frame, 
respectively. One obtains 

^ ^ 1 / iws UJ2 + iuJi\ ^ /a'S - P'-/ -a' (3 + (3'ot\ 



and 



2 + 1^1 J \-ia' + ai -7/3' + a5' ^ ' 



2.2 Mobius transformations 

There is another way to establish the SO (3) action of SU(2) via Mobius transformations 
of the Riemann sphere. It is summarized here for reference in Section For more detail 
see 1^, pp. 29ff]. 

The unit sphere S"^ = {(Cj^jC) G | + + = 1} is mapped conformally onto 
the extended complex plane C U {00} by stereographical projection from the north pole. 
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To the point z = x + iy in the complex plane corresponds the point in the sphere with 
coordinates 

2x 2y ^ 



>'-||2^T ^ = I 12^1 ' C= 2^-1 ■ (2.4) 

The conformal maps of the Riemann sphere onto itself are described by Mobius transfor- 
mations 

0-2: -)- p 

z , where a5 — = 1. (2-5) 

72 + 

The map sending a matrix in SL(2,C) with entries a, /3, 7, 5 to the Mobius transfor- 
mation ( p.5| ) is a group homomorphism with kernel {1,-1}. In particular, the isometric 
transformations of the sphere, i.e. the rotations, are conformal. They correspond to those 
Mobius transformations with 5 = a, 7 = — /3. This is the image of SU(2) under the above 
homomorphism. 



3 The Lagrangian, equations of motion, and their solution 
in terms of hyperelUptic functions 

The Lagrangian description of the toy top is used to derive equations of motion. The 
integrals of motion — total energy and two further integrals connected to the rotational 
symmetries of the system — are used to reduce the system to one degree of freedom. It 
is then solved in hyperelliptic integrals. Most of the results in this section, but not the 
equations of motion (^]^) and their derivation, are already contained in the classical works 
mentioned in the introduction. 



3.1 The Lagrangian, first integrals, and reduced equation of motion 

In general, the Lagrangian is a function on the tangent bundle of the configuration space. 
If this space is a Lie group, it is convenient to trivialize the tangent bundle via left or right 
multiplication. We will use right multiplication: 

SU(2) X su(2) ^ TSU(2) 
(<I>,u;) ^ (<I>,u;<I>). 

This corresponds to a description in the fixed frame. For a general exposition of this 
approach to mechanical systems similar to the Lagrange top, see Q, 

In these fixed frame coordinates, the kinetic and potential energy functions are 

r(c^, u:) = ^A{u;, io) + ^{C- A){u;, rf + ^ ps{[r, k],u;f, 

and 

=p{r,k), 

where /c = 63 is the unit vector pointing vertically upward, r = ^k^^^ is the unit vector 
pointing in the direction of the top's axis, s is the distance between the tip of the top and 
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its center of mass, p is the product of s and the mass, and C and A are the inertia moments 
of the top with respect to the symmetry axis and any perpendicular axis through the center 
of mass. The gravitational acceleration is assumed to be one, which can be achieved by 
a suitable choice of units, e.g. for mass. Note that r' = [io,r]. The kinetic energy of the 
toy top is composed of two parts: one (the first two terms in equation (|3.1D) is due to the 
rotation of the top around its center of mass and the other is due to the vertical motion of 
the center of mass. If this last term (ps/2)([r, k],uj)'^ was missing, and the inertia moments 
were related to the tip and not the center of mass, one would obtain the kinetic energy of 
Lagrange's top. The Lagrange function is 

L = T-V, 

with corresponding momentum 
dL 

m = — — = Alo + (C — A){lj, r)r + ps([r, A;], (^)[r, k]. 

OOJ 

It is convenient to introduce a derivative DL{(^^ijj) G su(2) with respect to the first 
argument of L by 



(L»L($,a;),r/) = ^L(e^''cI>,^) 



le=0' 



One obtains 

DL{^,uj) = {C - A){uj,r)[r,uj] +ps([r, k\,uj)[r, [k,uj]] +p[r, k]. 

Proposition 1. The equations of motion are 
m = lu), ml + DL(^, uj) 

(3.1) 

The first integrals 

h = T + V, l = {m,k), and n = {m,r) 
allow the reduction to 

where u = (r, k) . 

Proof. Consider a variation ^{t, e) of <I>(t, 0) with fixed end-points t = tQ, t = ti. Suppose 
d^/dt = Lo^ (yielding the second equation of motion) and d^/de = r/$. For the variation 
of the action functional S = J^^ Ldt one obtains. 



dS_ 
de 



ti 

{DL{^, uj) - [m, uj] - m', rj) dt. 
e=0 J to 
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(Note that duj/de — drj/dt + [uj,r]] = 0.) This yields the first equation of motion. The 
total energy h is always a first integral; / and n are easily checked by direct calculation. 
The integral I is due to the rotational symmetry around the vertical axis and n is due to 
the rotational symmetry around the top's axis. 

To achieve the reduction, note that in terms of the angular velocity, the momentum 
integrals are I = A{uj,k) + {C — A){u!,r){r, k) and n = C{uj,r). Thus one obtains {ui,r) = 
n/C and {uj,k) = 1/A + {1/A— 1/C)nu. Furthermore, {uj,[r,k]) = u' . Unless r = zt/c, 
one finds the following representation of to in the basis (r, k, [r,k]): 

(1 -u^)lo= (C-^(1 - u2)n - A-^{1 - nu)u) r + A-^{1 - nu)k + n'[r, k]. 

Now one can express the energy integral h in terms of /,n,ti and u' . This yields the reduced 
equation of motion (|3.2|) . ■ 



3.2 New dynamical constants 

We introduce new dynamical constants ei, 62, 63 to replace n, h and /. Together with 
±64, also introduced in this section, they are the branchpoints of the hyperelliptic Riemann 
surface defined in the next section. 



Proposition 2. The third degree polynomial on the right hand side of equation (3.1) has 



three real roots. Two of them lie in the closed interval [—1,1] and the third is greater or 



equal to 1 . Hence, equation (3.i) can he written 



_ ^ L.2 „2\ „./2 

where 



^yU^-el)u' = (u - ei) (n - 62) (li - 63) , (3.3) 



-1 < ei < 62 < 1 < 63 (3.4) 

and 




64 = W1 + — . (3.5) 



Proof. This follows from the fact that for the dynamical constants I, n and h to be 
physically feasible, there has to be a value for u between 1 and —1 such that equation 
(3.2) leads to real u' . This means that the polynomial on the right hand side must be 



nonnegative somewhere in the interval [—1,1]. But at u = ±1 it takes the nonpositive 
values —{I =F ^)^/(2^p), and it goes to itcxD for u — > ±00. Hence there have to be three 
real zeroes, situated as stated. ■ 



Conversely, it can be shown that inequality (3^) is the only constraint for the zeroes. 
I.e., given three numbers satisfying this inequality, there is a state of motion of the top 
with 61 < u < 62 that leads to exactly these zeroes. 

Note that the initial value problem with differential equation (^^) and an initial value 
for u between 61 and 62 has more than one solution. First, the initial sign of u' has to be 
determined. But even then the solution is only unique up to the moment when u reaches 
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Figure 1. The top sheet of the hypereUiptic curve C. 



ei or 62- After that, u can either reverse its path straight away, or stay constant for 
some time. Of course, this indeterminacy is not inherent in the original system, but is 
introduced by the reduction. Upon closer examination of the system, one finds that u is 



constantly equal to e only if e is a double zero of the right hand side of (3^). See [Q, pp 



280ff] for a discussion of the analogous situation in the case of the Lagrange top. Except 



for this singular case, which will be treated in Sections 5^ and p.3[ only those solutions 
of ( |3.3| ) have to be considered that oscillate periodically between the isolated minima and 
maxima ei and 62- 

From now on, we will use ei, 62 and 63 as dynamical constants instead of /, n and h. 
However, the first do not uniquely determine the latter. Indeed, substituting 1, —1 and 



into the right hand sides of (|3.2| ) and (^), one obtains the relations 
(l-ei)(l-e2)(l-e3) = -^(/-n)2 

(-1 - ei) (-1 - 62) (-1 - ea) = -^{l + nf (3.6) 
1 / n2 /2 X 

If the constants 61,62,63 are given and none of them is equal to 1 or —1, then the first 
two equations of have four solutions for (/, n). If (a, b) is one of them, the others are 
(—a, —6), (6, a) and (—6, —a). Note that replacing / and n by —I and — n is equivalent to 
looking at a mirror image of the system. Once a solution for (Z,n) is chosen, the third 
equation of ( [j.6| ) determines the value of h. Unless A = C, the four solutions for (/, n) 
lead to two different values for h. 

If one of the constants 61, 62, 63 is equal to 1 or —1, then I = n 01 I = — n, respectively, 
so that there are only two solutions. 

3.3 The hypereUiptic time integral 

Solving equation (|3.3|) by separation of variables, one obtains for t the hypereUiptic integral 



" ('^^ - dn. (3.7) 



2 (li - ei) {u - 62) {u - 63) 



Allow complex values for u and t so that it becomes an Abelian integral on a hypereUiptic 
Riemann surface. Let C be the hypereUiptic curve given by w'^ = R{u), where 
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R{u) = {u- ei){u - e2){u - e3)(u^ - el). 
For future reference we remark that by comparison with equations (^]^) and ( |3.3D , 

The hyperelhptic curve C is a two sheeted branched covering of the u-plane with the six 
branch points ei, 62, 63, 64, —64, and 00. Figure |l] shows one sheet of C with cuts along 
the real axis. This will be called the top sheet. Lift the path along which u moves during 
the motion of the top to a path on the hyperelhptic curve which is homologous to the path 
drawn in the figure. Then 




This is an Abelian integral of the second kind which has a simple pole at infinity: Introduc- 
ing a holomorphic parameter = 1/u around 00, one obtains the asymptotic expansion 
t = -iV2sv-^ + 0(1) for V ^0. 



3.4 Hyperelliptic integrals for the Cayley-Klein parameters 

In the following proposition, hyperelliptic integrals for the Cayley-Klein parameters are 
presented and their analytical properties discussed. These solutions are similar to the 
corresponding results for the Lagrange top. In the latter case, one obtains elliptic inte- 
grals and not hyperelliptic ones as here, but they have the same kind of singularities at 
corresponding places. Also, the reduction to the case ^ = C of spherical tops works for 
the Lagrange top as well. According to F. Klein p. 234], the reduction of Lagrange's 
top to the spherical case was first noticed by Darboux. 

Proposition 3. Suppose first that A = C . Then the solution of the system is given by 
the hyperelliptic integrals 

f 1 / 1*2-642 \ du 

J 2{U + 1) V 1-64^ J W 

J 2{U-1) \ 1-64^ J W 



n2 


-64^ 


1 


-64^ 






n2 


-64^ 


1 


-64^ 






n2 




1 


-642 






n2 




1 


-642 



, , 1 / m2 — 642 \ dn ^ ^ 

2[u — 1) V 1 — 64"^ J w 

1 / 1*2-642 \ du 

2{u + 1) V 1 — 64^ J w 

The constants w±i denote one of the two values of w on C over u = ±1, i.e. 



w+i = ±y/R{l) and W-i = ±y/R{-l) . 
(There are four possible ways to choose, in accordance with the indeterminacy of I and n 



in terms of ei , 62 , ea ; see Section 3.i) 
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Table 1. The poles of the logarithmic differentials. 



Differential 


pole with res. 1 at 


pole with res. -1 at 


daja 


(u,w) = (-1, 


oo 




(u,w) = {+l,-w+i) 


oo 


dlh 


{u,w) = (+1, W+l) 


oo 


dS/S 


{u,w) = (-1,-W-i) 


oo 



They are Abelian integrals of the third kind, the differentials under the integral sign 
having two simple poles each, with residues ±1 at the places shown in Table ||, and no 
other singularities. 

If A ^ C, the solution is 

where 

1 1 



C A 



n. 



I.e., the solution differs from the solution for the spherical top only by a rotation with 
constant speed around the top's axis. 

Proof. First, the expressions for the logarithmic differentials of a, l3,j,5 are derived in 
the general case. The reduction to the case of spherical tops is then immediate. Finally, 
the analytic properties of the differentials are examined. 
Observe that 

1 = aS — f3j ^ ^ 

(3.12) 

u = a6 + /?7. 
This implies 

u + l = 2a6 and u - 1 = 2/^7, (3.13) 

such that 

du = 2(5 da + a d5) 

^ ' 3.14 

du = 2{-fd(3 + pdj). 



By equations ( |2.2| ) and (^^) , the components of the angular velocity vector in the direction 
of the vertical axis and the top's symmetry axis are given by 

iuj3 = 2{a'6 - 13' -f) = 2{-a5' + f3-f'), 
iCls = 2{6a' - 13-i') = 2{-a5' + p'-f). 

This implies 

i{^z + oJz)dt = 2(5 da — a d5), 
i{n3-uj3)dt = 2{jdp- Pd'j). 
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With equations (|t|) and one obtains 



du + 


+ uj3)dt = 


A 6 da 


= 2{u + 


a 


du + i{i^3 


— uj'i)dt = 


A-fdf3 


= 2{u - 


dl3 
' P 


du — 2(1^3 


— uj3)dt = 


4/3 d7 


= 2{u - 


1)^ 

7 


du — 


+ uJ3)dt = 


4 a 


= 2{u + 


4, 



Now it fohows from n = CVt^ and / = Alo-^ + (C — A)Q.3U that 
n3+uJ3 = ^{l + n) + {u+l){^^- ^ n, 
- = -^{l -n)-{u-l)(^- ^ n. 
Note that from equation (^), 



R{±1) 



1 



2jrs 



Hence the signs of w±\ can be chosen such that 
V 2 



or, from equation (3.5), 



A 



Together with equation (3.16) and the equation for dt obtained from (| 

u — 64^ du N / 1 1 \ , 

?(i23 + io-Adt = ^ w-\ h z(tt + 1) -p: 7 \ndt 

w \L A 



-z(r23 — uj3)dt 



1-642 
u — 64^ 



w+i h t[u — 1) 

w 



1 1 

C~A 



Substituting these expressions into ( ^.15 ), one obtains 

du i ( \ 1 



da. 




^^w + 


u" 


-64^ 


a 


2(u + l) 


1 


-642 


dp 




w — 


u^ 




T 


2(u- 1) 


1 


-642 


dj 








-64^ 


7 


2(u- 1) 


1 


-642 


d6 




w — 




-64^ 


T 


2(u + l) 


1 


-642 



^ 2 ^ 
du ^ / 1 1 
F ~ 2 I C ~ I 



n dt. 

n dt 
n dt 



du i f ^ 1 , 
^-^^'F + 2U"I'"^' 

, ^ / 1 1 1 

"-^•f-21c-I'"^*- 



(3.15) 



(3.16) 
(3.17) 



this yields 
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This implies equations ( 3.10| ) and the reduction to the spherical case. 



Now assume that A = C. The assertions regarding the poles of the logarithmic differ- 
entials away from oo follow straightforwardly. Regarding the asymptotics at oo, note that 
near oo the logarithmic differentials are 

du , , , , 

+ a holomorpruc part. 



2(n±l) 

Introduce the parameter = ± 1). It is well defined (up to sign) and holomorphic 
around oo. One finds that the singular part of the logarithmic differentials \s d^/^. ■ 



4 The trajectory of the tip of the top 

In this section we examine the curve that the tip of the top describes in the horizontal 
plane. This aspect of the top's motion is particularly easy to observe experimentally. For 
example, in toy shops one can buy tops that have a pencil lead for a tip or a felt tip pen. 
Lord Kelvin did experiments with such pencil tipped tops, while F. Klein used cogwheels 
from a clockwork which he spun on a sooty glass plate. See pp. 619ff] for pictures and 
a discussion of the significant influence of friction and other concomitants neglected in the 
mathematical model. 



4.1 Geometrical considerations 

Since r is the unit vector in the direction of the top's axis, pointing away from the top's tip 
if attached to the center of gravity, the curve traced by the tip of the top on the supporting 
plane is the orthogonal projection of —sr onto that plane. The following proposition gives 
a formula for this curve in terms of the Cayley-Klein parameters. 

Proposition 4. Identify the supporting plane with the complex number plane. Then the 
curve traced by the tip of the top is 

c = 2sa(3. 

Proof. First, project the vector — r stereographically from the north pole of the unit 
sphere into the complex plane. The result is the image of under the Mobius transfor- 
mation ( |2.5| ), i.e. z = (5/5. Since z = — 7/a, the fromula for c follows from equations 
&■ 

2z 

c=s{^ + iri) = s 

p| + 1 

This proves the proposition. ■ 

For the qualitative considerations below, it is useful to represent the curve c in polar 
coordinates: c = pe^'^. Clearly, p = sVl — u'^- The angle cf) is not uniquely defined if p 
vanishes, i.e. if -u = ±1. Otherwise, the following proposition gives a formula for d(j). 



Proposition 5. If —1 < u < 1, then 

. 1 — 64^ / W-i w+i \ du 



2 1 — 64"^ V« + l u — 1 J w 



(4.1) 
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Proof. The angle </> is the argument of the curve z = 2sa(3. It follows that d(j) is the 
imaginary part of the logarithmic differential dc/c = da/a + d(3/p. Formula ( |4.1| ) then 
follows from equations ( p.lO ). ■ 



4.2 Loop, cusp, or wobbly arc? 

Figures ^ and |3| show different trajectories of the top's tip. Three qualitatively different 
cases are clearly discernible. In Figure ^ and the top row of Figure ^, the curve is a smooth 
line circling the origin. The first picture on the bottom row of Figure |3| shows a cusp, and 
the following two show loops. The next proposition gives a condition in terms of ei, 62 
and 63 for the three cases to occur. 

Proposition 6. Assume that —1 < ei and 62 < 1. If w^i and w^i have the same sign, 
then the trajectory of the top's tip has neither loops nor cusps. If Wj^i and W-i have 
different signs, then the curve has loops if and only if 

1 - 6162 
< 63. 

62 - ei 

There are cusps if both sides are equal. 

Proof. Loops occur, if d(p/du changes sign as u moves from ei to 62. It follows from 
equation ( [4.1D that this happens if W-i/{u + 1) — Wj^i/{u — 1) = has a solution for u in 
the interval (61,62). That equation is equivalent to 

^ = ^ . (4.2) 

U — 1 Wj^i 

Since the left hand side of this equation is smaller than zero for u S (—1,1), this equation 
cannot be fulfilled if Wj^i and w^i have the same sign. This proves the first part of the 
proposition. 

Now assume that Wj^i and w_i have different signs. Then 



w-i _ _ / (-1 - 6i)(-l - 62)(-l - 63) 

W+l~ V (l-6i)(l-62)(l-63) 
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Figure 3. The curve traced by the tip of the top. Here, W-i and w+i have different signs, 
s^l, ei = 0.7, 62 = 0.9, = 2 and 63 = 1, 1.1, 1.3, 1.85, 2.5, 3 (top left to bottom right). 

In the interval [61,62], the function u ^ {u + l)/(u — 1) is continuously decreasing from 
(61 + l)/(ei — 1) to (e2 + l)/(62 — !)• Hence, equation ( |4.2| ) leads to a n G (61,62) if 

62 + 1 _ / (-l-ei)(-l-e2)(-l-e"^ 61 + 1 
62-1 V (l-ei)(l-e2)(l-63) 61 -r 

This is equivalent to 

1 - 63 ^ 1 + 61 1 - 62 ^ 1 + 63 
1 + 63 1 - 61 1 + 62 1 - 63 ' 

The inequality on the right is always fulfilled, because the right hand side is greater 
than one and the left hand side smaller. The inequality on the left is equivalent to (1 — 
6162)7(62 — 61) < 63. This proves the condition for loops. If (1 — 6162)7(62 — 61) = 63, 
then the zero of d(j)/du occurs at u = 62, hence there are cusps. ■ 



5 The degenerate cases 

If two zeroes of the polynomial R{u) coincide, the hyperelliptic integrals of the Sections 



3.3 and 3.4 degenerate to elliptic ones. In those cases, it is possible to solve the system in 
terms of elliptic functions. This will be done in the following sections. Since — 1 < 61 < 
62 < 1 < 63 and 64 > 1, only the following cases have to be considered: ei = 62, 62 = 63 
and 63 = 64. Furthermore, the limit case 64 = —64 = 00 is examined, because it is of 
special interest: In this case the toy top becomes the Lagrange top. 
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5.1 The case 63 = 64 

Suppose that 63 = 64. Then the hyperelhptic curve C degenerates to the ehiptic curve 
Ce3=e4, given by = R{u), where 

R{u) = 4(n - ei)(n - e2){u + 63) . 



The leading coefficient of R is chosen to be 4 in accordance with the Weierstrass normal- 
ization of elliptic curves. This will simplify the integration of the elliptic integrals in terms 
of Weierstass elliptic functions p, ( and a. The branchpoints are ei, 62, —63, and 00. 
The i-integral and the integrals ( 3.10| ) become 

t = -iV2^ [ (u + es)— , (5.1) 



and 



, 1 / M + 63 \ du 
^°Sa= / — — — w + —— w-i — 

2{U + 1) V -1 + 63 J W 

1 / n + e3 \ du 
logp = / — — w - — w+i 



2{u-l) V 1 + 63 "V (5 2) 

f 1 / u + e^ \du 

J 2[u-l) \ 1 + 63 J w 

f 1 / u + es \ du 

J 2{u + l) \ -1 + 63 J w 

As before, w±i denotes one of the two values of w on the curve Ce3=e4 over u = ±1, i.e 



U7+1 = ±a/R(1) and w^i = ±^/R{-l). 
5.1.1 Uniformization of the eUiptic curve 

The system will be solved by using an Abelian integral of the first kind to uniformize the 
the elliptic curve Ce3=e4 and then expressing t and a, j3, 7, 5 as functions of the uniformizing 
variable. 

The left part of Figure |^ shows the top sheet of the curve Ce^=ei with cuts and a normal 
cycle basis a, (3. Let wi and 002 be the corresponding half periods: 

f du f du 

2wi = / — , 2u2= —. 

Ja W Jp W 

Note that oji lies on the negative imaginary axis and that 0J2 is positive. Let the uni- 
formizing variable be given by 

ru=u(x) J 

x= — . (5.3) 

Ju=e2 W 

One obtains u and w as doubly periodic functions of x with periods 2loi and 2^2, namely 
n(x) = p(x - W2) + ^(61 + 62 - 63) (5.4) 
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Figure 4. Left: The top sheet of the eUiptic curve 6*63=64 • Right: A fundamental domain in the 
a;-plane. 



and 

w{x) = p'{x - UJ2). (5.5) 



The additive constant (ei+e2— e3)/3 in ( |5.4D comes from the fact that the branchpoints are 
not centered around zero as required by the Weierstrass normahzation, and the p-function 
is shifted by wi because the integral in ( ^ ) does not start in 00 but 62- 

The right part of Figure |^ shows a fundamental rectangle in the x-plane. The points 
X = iba and x = ±6 will be of importance in Section 5.1.3| and are defined as follows: 



The point a in the x-plane is supposed to correspond to the point {u,w) = {—l,W-i) on 
Ce3=e4, and b is supposed to correspond to (n, w) = (+1, tw+i). For example, let 



and 



a = ±^uji + j \JR{u) diL^ 
b = ± f ^R{u) du, 

J eo 



where the signs are chosen according to whether the points {u,w±i) lie in the bottom or 
top sheet. (The figure shows the case where both x = a and x = h correspond to points 
in the bottom sheet. Otherwise exchange +a with —a, or +6 with —6, respectively.) 

Regarding the toy top we adopt the convention that during its motion, the correspond- 
ing point on Ce3=e4 moves on a path homotopic to the cycle a. (This choice determines 
the sign on the right hand side of equation ( ^.ID ). The corresponding point in the x-plane 
then moves on the imaginary axis in the negative direction. 



5.1.2 Solution for time as function of the uniformizing variable 

We will now pull back the t-integral (|5.l| ) to the x plane and solve it. The following 
proposition gives the resulting formula for t as function of x. The initial condition which 
is chosen means that n = 62 at t = 0, i.e. the axis of the top is initially in its most upright 
position. 
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Proposition 7. Assume that t = for x = 0. Then equation {5A) implies 



t = i\/2s{C,{x - UJ2) + C{'^2) - -^(ei + 62 + 263) x) . 



Proof. Equations (^jj) and ( [5.5D imply that the pullback of the holomorphic differential 
du/w to the x-plane is dx. The t- integral therefore becomes 

t = -i^Pls J - UJ2) + ^(ei + 62 + 2e3))dx. 
Since C' = ~P; this implies 

t = iV2s[C,{x — W2) — -(61 + 62 + 2e3)x + const.). 
The constant has to be C('^2) for t to vanish at x = 0. I 

5.1.3 Solution for the Cayley-Klein parameters as functions of the uniformiz- 
ing variable 

The following proposition gives formulas for a, (3, 7, 5 as functions of x. But first, a few 
words have to be said about the initial conditions oq, Po, 10,^0 for t = a; = 0. 

It is no essential restriction to assume that ao is real and nonnegative and therefore 
equals 60, and that 70 is purely imaginary with nonnegative imaginary part, such that 

= 7o- For one may always achieve this by suitable rotations of the system about the 
z-Axis and the top's symmetry axis. 

Proposition 8. Under the above initial conditions, one obtains 
I a{x - a) 



a[x - UJ2) 

a{x + L02) 

7 = ^3e'3^^^, 
a{x - UJ2) 

a{x + 0)2} 



(5.6) 



where 



, , ,1 + 62 cf{^2) 

ki = ki = \ — — , 

2 a[a) 



, , ■ / I - 62 c^(^2) 



(5.7) 



and 



W—l 

h = -k = TT^-r— — r + C(a - ^^2), 

h = -^^ = ^(^+<(b-U.2). 
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Proof. Consider the asymptotic behavior of the logarithmic differential of a: 

da 1 ( It + 63 \ du 
w H W-\ — . 



a 2(u + l) \ —1 + 63 ) w 

It has only two simple poles on Ce3=e4, one at (n, vS) = (—1, w^\) and the other at 00. For 
(u^w) — > {—l,W-i), the asymptotic expansion is 

da f I , 
+ 0(1) du. 



a \u + 1 

To calculate the asymptotic expansion around the branchpoint 00, introduce the local 
parameter ^■^ = 1/u. The result is 

Since the residues of the logarithmic differential of a are 1 and —1, respectively, a itself is 
not branched on 6*63=64, but has a zero at (n, = {—l,w-i) and a simple pole at 00. It 
is a multiply valued function, because the additive periods of the logarithmic differential 
lead to multiplicative periods of a. It follows that, in terms of the uniformizing variable 
X, a is of the form 

a = fcie^ — ^ -. 

aix - UJ2) 



Similarly, one obtains the other equations (5.6). 



The constants fci, . . . , ^4 and Zi, . . . , Z4 are determined by the initial conditions. Note 
first that ao = Sq implies ki = and Pq = 70 implies k2 = k^. Just substitute x = into 
equations ( |5.6D and observe that a is an odd function. 

Further, since u = 62 for a; = 0, it follows from equations ( |3.13| ) that 

62 + 1 = 2ao5o = 2k^ ^'^"^ 



0-2 (a;2) 
and 



62 - 1 = 2/3o7o = 2h 



2 

From this one obtains equations ( |5.7D . The signs of ki and k2 are determined so that a^ 
is positive and 70 has a positive imaginary part. 
Concerning the constants h, ■ ■ ■ ,U, note first that 

^ = a6 = k^k,e(^^+^^^^ a(x - aMx + a) 



a{x - 0J2)cr{x + 002) 



is a doubly periodic function of x. This implies /4 = —li since the quotient of cr-functions 
is already doubly periodic. Equally, considering (3'y = {u — l)/2 yields I2 = —l-s- 
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From the logarithmic derivative of a with respect to x 



dloga a'{x — a) a'{x — ijJ2) 

~ ti + 



dx 

one obtains 

d\og a 



u 



dx 



a{x — a) a{x — UJ2) 



({x - a) + ({x - UJ2), 



since a' /a = The first of equations then ( |5.8| ) follows from 



lim 



dlog a 



C(x - a) 



dx ' 7 2(-l + e3) 
Since C{x — a) = l/{x — a) + 0{x — a) as x ^ a, one has to show that 



dlog a 



1 



+ 



W-l 



+ 0{x - a) 



dx X — a 2(— 1 + 63) 
From the first equation (p.2|) and dx = du/w one obtains 



(5.9) 



dlog a 



1 



n + 63 
w H : W-i . 



dx 2{u+l) y ' -1 + 63 
Since w = du/dx, this yields the expansion 



dlog a 

dx 



1 



du 



2(u + 1) V dx 
1 du 



+ {x - a) 



fu 
dx^ 



+ 0(x-aY+{l + 



u + 1 

-1 + 63 



u + I dx 



+ 



dx'' 



2 — 

x=a dx 



+ "f"' . + Oix - a) 



2(-l + e3) 



(5.10) 



The last equality uses 
X — a 1 



■u+1 

dx I 



+ 0{x - a). 



(5.11) 



The asymptotical expansion (|]9D follows from ( p^ and 



1 



1 du 



x — a u + \ dx 



1 

+ - 



d2n 
dx'^ 



2 — 

^=0- dx 



+ 0{x - a). 



To see this last equation, just take the Taylor expansion 

(x - af + 0{x - af, 



du 

u + 1 = — 

dx 



. 1 d^u 
(^-") + 2 d^ 



divide by {x — a) and (u + 1), and use ( |5.1l| ) again. The equation for /3 is derived 
analogously. ■ 
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5.2 The case of regular precession 

If a spinning top moves in such a way that the angle between the top's axis and the vertical 
is constant, its motion is called regular precession. It then follows from the symmetries of 
the system that the top spins around its axis with constant speed, while the axis precedes 
uniformly around the vertical axis. Notably, the solutions for t and the Cayley-Klein 
parameters in terms of hyperelliptic integrals fail in the case of regular precession because 
the path of integration degenerates to a point. One would have to look at the system 
before it is completely reduced to analyze this case. It turns out that regular precession 
occurs if R{u) has a double zero at the constant value for u. There are two qualitatively 
different ways this can happen. If ei = 62, small perturbations of the system will lead to 
small nutation. The regular precession is stable. But if 62 = 63 = 1, small perturbations 
will lead to nutation between u = ei and u = 62- The regular precession is unstable. See 
10, pp. 278ff] for an analogous discussion of the stability of regular precession in the case 
of Lagrange's top. 



5.3 The aperiodic case: 62 = 63 = 1 

If 62 = 63, there are two possible types of motion. Either u is constantly 1. This is the 
unstable case of regular precession mentioned above. Or u moves from 1 to ei and back. 
This case will be considered in this section. We proceed similarly as in Section ^.Ij . 

The hyperelliptic curve C degenerates to the elliptic curve Caper given by ■w'^ = R{u) 
with 

R{u) = 4(u - ei)(u2 - 64^). 
Branchpoints are ei, ±64 and 00. The t-integral becomes 

t = ^V2-s f^^l^^. (5.12) 
J u — 1 w 

It is singular at u = 1. That is why this is the aperiodic case. The integrals ( |3.10| ) become 

f 1 / 2 m2 _ g^2 \ 

log a = / — w ttW^i — 

^ y 2(u + l) V u-1 1-64^ J w 



log/? 
log 7 



du 



2(u- 1) 

du 
2{u-l) 



(5.13) 



, , , 1 / 2 n2 - 64^ \ du 
logo = / — w H ttW-i — . 

^ ' 2{U+1)\ U - 1 1 - 64^ J W 



5.3.1 Uniformization of the elliptic curve 

Figure ^ shows one sheet of Caper, which will be called the top sheet. Proceed as in Section 
|5Xl| , but let the curves c start in (u, w) = (ei, 0). For a normalized cycle basis, take a path 
a in the top sheet which encircles ei and 64 in the counterclockwise direction, and choose 
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-64 




w = - i I R 



1/2 I 



w = + 



1 


u=l 

X=C02 
u=-e4 


• 1 

x-a 

u=-l 

1 


x--a 
u=-l 

,x=b 
u=l 



X-O) 



x=(03 
u=oo 



Figure 5. Left: The top sheet of the elhptic curve Caper- Right: A fundamental domain in the 
a;-plane. 



P accordingly. Let uji and 102 be the corresponding half periods and let ^3 = uji+uj2- One 
obtains the uniformization 

U = p{x- UJ3) + y 

and w = p'{x — W3). 

The right part of Figure |5| shows a fundamental domain in the x-plane. Again, the 
points above u = ±1 will be important. In the figure, they are marked by asterisks. Let a 
and h be points in the x-plane corresponding to {u, w) = (—1, W-i) and {u, w) = (+1, w+i), 
for example. 



\/ R{u) du and b = ±i J \/—R{u) du. 



(In the figure they are drawn for the case in which (—1, W-i) lies in the top sheet and Wj^i 
has positive imaginary part.) 

Let the path of integration corresponding to the motion of the top be homotopic the 
the path with arrows drawn in the figure. 

5.3.2 Solution for time as function of the uniformizing variable 

The integrand of the t-integral has simple poles over n = 1 which leads to logarithmic type 
singularities for t. The Riemann surface of the function t will therefore be an infinitely 
sheeted cover of Caper with the two branchpoints (n, w) = (1, itw+i). Place a cut on Caper 
along the dotted line in Figure ^ and cut the x-plane along corresponding lines. Then 
our path of integration on Caper does not cross the cut, so that we can consider only one 
branch of the function t on the cut elliptic curve, or the cut x-plane, respectively. 

Proposition 9. Suppose that t = at x = {). Then 
* C{x-uj^)-C{oo^) + (l + ^-)x 



i^/2s ' ' ' \ ?, 
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where that branch of the function 

is chosen that is zero for x = 0. 
Proof. Put equation ( 5.12| ) in the form 

t = ^V2-s [(u + l+'-^)^ (5.14) 

and consider separately the integrals 
1 du 



h 
and 



u — 1 w 



/, du 
{u + l) — . 
w 

Substituting the uniformizing variable in the first integral, we get 

dx 

h 



1 f p'{b — ui^) dx 



p'{b-uJ3)J p(x - W3) - p(6 - Us) ■ 
Using the formula 

p'iv) 



p(0 - piv) 

one obtains 



C(e-r/)-C(C + r?)+2C(r/), 



h = r I (C(x -b)-a^ + b- 2uj^) + 2C(6 - u;^))dx 

= — i / (C(x-6)-C(x + 6) + 2C(a;3) + 2C(6-^3))d2;. 

p'[b - W3) y 

The last equality follows from ("(^ + 2uJk) = C(C) + '^Ci^kj- Now the last integral can easily 
be solved since = ^'iO/'^iO- Note that 

- a;3) = w+i 

to obtain 

h = ^ (log ( ) + 2(C(a;3) + C{b - ^a))^ + const) . 

w+i V \a{x + b)J J 

Requiring that /i = for x = 0, one obtains 
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where that branch of the logarithmic term is chosen that vanishes for x = 0. The other 
integral is easier to deal with: 

l2 = j (p(x-W3) + l + ei/3) 

= —C{x — UJ3) + (1 + ei/3)x + const., 
since p = — C'. If the constant is chosen to make I2 vanish for x = 0, this becomes 

I2 = -C(x - u;3) - C(W3) + f 1 + ^ V. (5.16) 



3 

The formula for t follows from equations ( 5.14[ ), ( |5.15| ) and ( ^.16 ). 



5.3.3 Solution for the Cayley-Klein parameters 

The functions a,/3, 7, (5 are also branched at the points above u = I. So we apply the 
same cut as in the previous section and have to choose branches. As in Section |5.1| we will 
assume that the initial conditions oq and 60 are real (and therefore equal) and positive, 
and that /?o and 70 are imaginary (and therefore equal) with positive imaginary part. 

Proposition 10. // the initial conditions are chosen as explained above, then 
^ _ i.>-i^P ^(^ - °) f(^ix + b)y 



a{x — uj^) \a{x — b) 



7 = /? = y 

^ ^ j^^-i^+i^p_(T{x + a) fa{x-b)Y 



where 

P 
k 



a{x + LJ3) \cr{x + b) 

1 W-l 

2 w+i ' 



1 + ei <y{i^-i) 
2 cr(a) '' 



4 1 — 64^ 

and those branches of the multiply valued functions {cr{x + b)/a{x — b)Y and 
((t(x — b)/a{x + b)Y are chosen which take the values e™^ and e"*'^^ at x = 

Proof. The formula for (3 and 7 follows elementarily from the corresponding integrals 
( 5.13| ) and the initial conditions. Now consider the logarithmic differential of a from 



equations ( |5.13| ): 

da f 1 1 w^i 1 W-i 



a \2{u + l) (ii-l)(u + l) w l-e^^ w 



du. 



As one can read off from this expression, the differential has three simple poles away from 
infinity, one at {u,w) = {—l,W-i) with residue 1 and two more at {u,w) = (+1, ittD+i) 
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Table 2. The simple poles of da/a and d5/5 with residues. 

da/ a 



{u,w) ^ 




00 




(+1, -w+i) 


Residue; 


1 


-1 


-w_i/2it;+i 


w_i/2w+i 






dS/S 






{u,w) = 


(-1,-^-1) 


00 






Residue: 


1 


-1 


w-i/2w+i 





with residues ^W-i/2w^i. Furthermore, the term du/2{u + 1) contributes one more 
simple pole at infinity with residue —1. The location of the poles with their residues is 
summarized in Table ^. It also lists the poles of the logarithmic differential of /3, which 
are obtained similarly. Now the poles with residue ±1 lead to zeroes and poles of a and 
6, while the other poles give rise to branchpoints. It follows that, as functions of x, a and 
6 are of the form 

p 



a 



P 



h e 



hx (^{x - a) [ (j{x + b) 
a{x — b) 
a{x — b) 



a[x - W3J 
a(x + a) 



a{x + UJ3) \a{x + b) 



(5.17) 



Choose the branches of these multiply valued functions as explained in the proposition. 
Then the values 



1 + ei a{uJ3) 



cr a 



1 + ei ajuJs) 
2 a{a) 



-iirp 



follow from the initial condition ao = Sq = {I + ei)/2. 

Since 2a6 = n + 1 is a doubly periodic function of x, it follows that I2 = —h- It is left 
to show that /i = / as given in the proposition. Since from equation ( ^.17 ) 

d log a 



dx 



h + C{x - a) - C{x - 003) +p{C{x + b)-C{x-b)), 



this will be achieved if we can prove that for x 



a, 



dlog a 



1 



+ 



W-l 



W-l 



+ 0{x-l) 



dx X — a 4 1 — 64^ 
This follows by a calculation analogous to the one in Section ■ 

5.4 The Lagrange top as limit case: 64 = —64 = 00 

Look at the t-integral as it is written in equation (3.7). If one simply lets 64 tend to 
infinity, the integrand diverges. But note that from equation ( |3.5| ) one obtains 

s _ A 1 

~2~ 2^ 1 - 64^ ' 
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so that equation (|3.7[) can be rewritten as 

du 




- ei){u - e2){u - 63) 



Now if one lets 64 tend to infinity while keeping A/p constant, this converges to the Abelian 
integral of the first kind 




on the elliptic curve w"^ = {u — ei){u — e2){u — 63). 

The integrals ( |3.10| ) cause no problems. They become 

f w + W-i du 
^oga = J , — 

log/3 = I 

log7 = j 

log6 = J 

This is exactly the solution F. Klein obtains for the Lagrange top p. 238], pp. 28f]. 
This is not surprising since this limiting case is obtained by letting s tend to zero while 
keeping all other parameters constant. But in this case the Lagrangian ( |3.lD coincides 
with the one of Lagrange's top. 



2(n+l) 


w 


w — W+l 


du 


2(n- 1) 


w 


W + W-^l 


du 


2(u- 1) 


w 


W — W-l 


du 


2(m + 1) 


w 
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